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ABSTRACT 


The study of mutualistic plant and animal networks is 
an emerging field of ecological research. We 


reviewed progress in this field over the past 30 years. 


While earlier studies mostly focused on network 
structure, stability, and biodiversity maintenance, 
recent studies have investigated the conservation 
implications of mutualistic networks, specifically the 
influence of invasive species and how networks 
respond to habitat loss. Current research has also 
focused on evolutionary questions including 
phylogenetic signal in networks, impact of networks 
on the coevolution of interacting partners, and 
network influences on the evolution of interacting 
species. We outline some directions for future 
research, particularly the evolution of specialization 
in mutualistic networks, and provide concrete 
recommendations for environmental managers. 
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INTRODUCTION 


Plant-animal mutualistic interactions, such as pollination and 
seed dispersal, have been regarded as one of the principle 
examples of coevolution and have long drawn the sustained 
attention of researchers (Fleming & Kress, 2013). The study of 
mutualistic interactions is part of the larger field of ecological 
networks, which also includes ‘antagonistic’ interactions 
between different species at varying trophic levels in the food 
chain (Pimm, 1982; Thompson, 2009). However, since the 
evolutionary forces that shape antagonistic networks may differ 
from those that shape mutualistic ones, we focused on 
mutualistic plant-animal interactions, a subject with now 
sufficient studies to stand on its own (Thompson, 2005; 
Thébault & Fontaine, 2010). Indeed, the published literature is 
so large and complex that is it useful to categorize the different 
kinds of questions that have been asked. Our review had the 
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following objectives: (1) outline the major ecological and 
evolutionary questions about mutualistic networks that have 
become prominent in the last 30 years, and (2) point out lines of 
research where further development would be particularly 
fruitful. ' 


BRIEF HISTORY OF RESEARCH ON PLANT-ANIMAL 
MUTUALISTIC NETWORKS 


The idea of plant-animal interaction networks was first proposed 
150 years ago with the publication of Darwin’s On the Origin of 
Species, who wrote “Il am tempted to give one more instance 
showing how plants and animals, most remote in the scale of 
nature, are bound together by a web of complex relations” 
(Darwin, 1859, p74). A large number of studies have looked at 
specific mutualistic interactions between plants and animals, 
such as yucca moths and yuccas (Pellmyr, 2003), ants and 
Acacia (Janzen, 1966) and fig wasps and figs (Cook & Rasplus, 
2003). However, the interactions among organisms are far 
more complex than pair-wise relationships, as many different 
organisms interact together in networks that include many 
species (Thompson, 2005). With more biological data available 
and the advance of algorithms for complex network analysis 
derived from physics and computer science, ecologists and 
evolutionary biologists have started to explore plant-animal 
mutualistic interactions from the network perspective (Guimaraes 
et al, 2011; Olesen et al, 2007; Rezende et al, 2007a, b). 

To study the growth in this literature, we searched the topics 
‘mutualistic network’, ‘plant animal interaction’ or ‘mutualistic 
interaction’ with the timespan from 1970 to 2013 on Web of 
Science, and then further refined the results using the Web of 
Science ecology category. A total of 1 945 publications were 
closely related to the topics, with rapid growth year-to-year 
(Figure 1). In terms of the total number of publications in a year, 
most years had more publications than the year before. This 
trend is expected to keep on increasing annually in the future. 
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Figure 1 Literature found on Web of Science by searching ‘mu- 
tualistic network’, ‘plant animal interaction’ or ‘mutualistic inter- 
action’ with the time span 1970 to 2013, plotted against the date 
of publication 

Each bar represents the number of publications published in a single 
year. This search was conducted in June 2014. 


Tracing early literature, plant-animal mutualistic studies 
began to appear in the late 1970s, with an escalation in 


publications observed in the early 1990s. Most early studies 
were field experiments, which laid the foundation for later, more 
theoretical research. Jordano (1987) detailed the first well- 
recognized study on mutualistic networks from the theoretical 
perspective in American Naturalist. A rapid increase in 
publications on networks was observed from 2005 (Figure 1). 
Bascompte et al (2003) introduced an important technique for 
network analysis to measure the nestedness of the interaction 
network, opening new avenues in the field. Beyond studies on 
network architecture, increasing research has focused on the 
implications of network theory for ecology and evolution, and 
we address these two maior kinds of studies separately below. 


MUTUALISTIC NETWORKS IN ECOLOGICAL TIME 


Mutualistic network properties 

Plant-animal mutualistic networks can be described by 
interaction matrices, with plant species in the columns and 
animal species in the rows; we refer to an interaction between a 
plant and animal species as a “link”. According to the kind of 
values in the matrix cells, network matrices can be categorized 
into “weighted” networks (“quantitative’ networks) and 
“unweighted” networks (“binary’ networks or “qualitative” 
networks, as they only indicate whether a pair of species 
interact, not the intensity). Several network properties have 
been generated to describe the network structure, such as 
connectance and interaction strength (Table 1). 





Measures the imbalance in the interaction strength of two interacting species (Bascompte et al, 2006). It is defined as 
ASj=(b;-bj)/(by+bj), where bj is the reciprocal dependence of species j on species i (bj; see interaction strength, Bas- 


Proportion of the realized interactions in all possible interactions (Yodzism 1980). In mutualistic networks, connectance (C) 
is: C=L/(IJ). L describes the number of realized links; / and J are the number of species of each bipartite network (Bluth- 


Interaction strength of species j on species i (bj) can be defined by the proportion of interactions between i and j (aj) of the 
total interactions recorded for i, thus p — DE. a,- For mutualistic network, b; measures the dependence of species i 
j ~ Aij ja Ý 


Measures the degree to which the network is organized into clearly delimited modules (Olesen et al, 2007). Modularity 
, where Nm is number of modules in the network, ls is the number of links in the network, and ks is 


Table 1 Glossary of main network terms and key references 

Metric Definition 
Asymmetry 

compte et al, 2006; Blüthgen et al, 2007). 
Binary network In binary network matrix, the value is 0 or 1, if the interaction occurs, the value is 1, otherwise 0 (Jordano, 1987). 
Connectance 

gen et al, 2008). 
Degree Number of interactions a species has (Jordano et al, 2003). 
Interaction 
strength 

on its partner j (Jordano, 1987; Blüthgen et al, 2007). 
Modularity 

(M O WI ko 

M= os (sy 
2 EGP 

the sum of degrees of all species in s. M values belong to the interval [0; -1/Ny]. 
Module A set of weakly interlinked subsets of species that consist of strongly connected species (Olesen et al, 2007). 
Nestedness 


Weighted network 


A nonrandom pattern of the network structure, which entails the tendency of specialized species to interact with a subset 
of the interaction partners of more generalized species. The nestedness temperature (7) measures the departure from a 
perfectly nested interaction matrix, ranging from 0 to 100, which indicates the degree of disorder. T=0 is defined for maxi- 
mum nestedness: when rows and columns are ordered by decreasing number of links, links of each row and column 
exactly represent a subset of the previous ones. Nestedness can be defined as N=(100°-7)/100° (Bascompte et al, 2003). 


Networks that include information on the intensity or weight of the interactions among nodes (Bascompte & Jordano, 
2007). 
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The topology of ecological interaction webs (such as 
mutualistic networks) holds important information for biodiversity, 
ecosystem stability and theories of coevolution (Bascompte et 
al, 2006; Jordano, 1987; Montoya et al, 2006). On the whole, 
mutualistic networks are neither a collection of pair-wise, highly 
specific interactions nor diffuse, random assemblages 
(Bascompte & Jordano, 2007). Rather, mutualistic networks 
have common, well-defined network architecture regardless of 
the type of mutualism, species composition, and geographic 
region (Bascompte et al, 2003). Several topological features are 
characteristic, including skewed distribution of links per species 
(i.e., a few species with many more interactions than expected 
by chance, and many species with a few interactions; Vazquez 
& Aizen, 2003; Jordano et al, 2003), the nested organization of 
the interaction matrix (Bascompte et al, 2003) and the frequent 
occurrence of asymmetric interactions (i.e., a plant species 
depending strongly on an animal species, the animal depending 
weakly on the plant; Bascompte et al, 2006). In addition to 
these properties, some networks are also modular (especially 
pollination networks; Olesen et al, 2007), whereby clusters of 
species interact more closely with each other than with species 
in other clusters or species outside the clusters. In comparison 
to antagonistic networks, mutualistic networks demonstrate high 
asymmetry and connectance, whereas antagonistic networks 
tend to be weakly connected, with many groups of species 
(“compartments”) interacting only within their group (Cagnolo et 
al, 2011; Krause et al, 2003; Prado & Lewinsohn, 2004; 
Thébault & Fontaine, 2010). 


Stability and diversity of mutualistic networks 

The maintenance of stability in complex communities has been 
a long-standing debate since the classic work of May (1974). 
Most network studies have found that the properties of 
mutualistic networks contribute to their diversity and stability; for 
example, the asymmetric nature (i.e., asymmetric interactions 
and asymmetry in interaction strength) of mutualistic networks 
promotes community coexistence, which favors the persistence 
of biodiversity (Bascompte et al, 2006). Nestedness has been 
shown to reduce interspecific competition, and hence promote 
diversity and nested networks stability (Bastolla et al, 2009; 
Thébault & Fontaine, 2010). Other network properties, including 
community size, species degree, species strength, and 
symmetry of the interaction, may also positively contribute to 
stability (Okuyama & Holland, 2008), while modularity 
decreases stability (Thébault & Fontaine, 2010). Recent 
developments have questioned the idea that nestedness itself 
affects stability, suggesting that more simple features, such as 
species degree, are more important drivers (Feng & Takemoto, 
2014; James et al, 2012; Jonhson et al, 2013; Staniczenko et al, 
2013). A further complication is that other interaction types, 
such as antagonistic networks, may interact with mutualistic 
networks to affect overall community stability (Mougi & Kondoh, 
2012). 


Spatial-temporal variation of the mutualistic network 
Studies have indicated that both spatial and temporal 
dimensions impose constraints on plant-animal mutualistic 


interactions and influence network patterns (Burkle & Alarcon, 
2011). Nevertheless, some general network properties remain 
constant over time and space. For example, Plein et al (2013) 
found that none of the characteristics of the seed-dispersal 
networks (e.g., interaction diversity, interaction evenness, and 
network specialization) used in their study changed with 
landscape type (farmland, orchard, forest edge) or season. 
Other studies have shown that although seasonal species 
turnover exists in the network, general network patterns (such 
as nestedness, connectance and modularity) remain relatively 
constant (Dupont et al, 2009; Dupont & Olesen, 2012; Olesen 
et al, 2008, 2011; Petanidou et al, 2008; Plein et al, 2013). 
However, these studies were conducted over a relatively short 
time period (often two to four years), and further studies with 
longer time spans are needed to test these conclusions. 
Additional attention needs to be paid when comparing different 
networks or when pooling data from different networks for meta- 
analysis. Networks are only comparable in this manner if they 
are at the same time scales or over the same phenological 
periods; for example, mixing networks from different seasons 
may prove confusing (Burkle & Alarcon, 2011). 


Conservation implications of mutualistic networks 

One potential problem for networks in our changing world is that 
alien mutualists, including both plants and pollinators, can 
integrate into native pollination networks, and sometimes end 
up acting as super-generalist species (a few species that 
interact with an extremely large number of species) of the 
network (Aizen et al, 2008; Bartomeus et al, 2008; Olesen et al, 
2002). Highly invaded networks exhibit weaker mutualism than 
less invaded networks, and the connectivity among native 
mutualists declines, although overall network connectivity may 
not change (Aizen et al, 2008) and other aspects of network 
structure such as nestedness are relatively robust to the 
introduction of invasive species (Vila et al, 2009). However, the 
removal of invasive alien species from the network can change 
the network structure, imposing a pronounced effect on degree 
distribution and modularity of the network, leading to higher 
species loss, which could affect the evolution of the interaction 
network architecture (Valdovinos et al, 2009). 

Habitat loss and fragmentation also have large, consistently 
negative effects on biodiversity (Fahrig, 2003). The extinction of 
species and loss of mutualistic partners may impose a 
cascading effect on mutualistic networks, which can 
consequently lead to greater biodiversity loss (Anderson et al, 
2011). Several theoretical studies have looked at how 
mutualistic networks respond to such disturbances. For 
example, Fortuna & Bascompte (2006) used a modeling 
method to understand how mutualistic networks respond to 
habitat loss, and found that real communities started to decay 
sooner than random communities, although they persisted 
better at high levels of destruction. Pollination network field 
research also demonstrated that habitat loss not only leads to 
species loss, but also indirectly causes the reorganization of 
interspecific interactions in the local community (Spiesman & 
Inouye, 2013). The reduction in suitable habitats is associated 
with species loss, which is correlated with reduced nestedness 
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and increased modularity (Tilman et al, 1994). 

Given the work on how mutualistic networks respond to 
invasive species and anthropogenic disturbance, it is clear that 
network theory has important implications for conservation. For 
example, by looking at the network structure, we can identify 
the network’s susceptibility to alien species invasion (Olesen et 
al, 2002; Morales & Aizen, 2006). Network theory can also 
clarify the role of different species (i.e., species with different 
specialization degree) and their impacts on the whole network 
architecture and stability. For example, super-generalist species 
can greatly affect the overall topology of a network (Aizen et al, 
2008; Hansen & Galetti, 2009), which has a direct consequence 
on the protection of the endangered system (Kiers et al, 2010) 
as these super-generalists can be targeted in conservation 
plans. Other issues that lead to species loss in the network, 
such as the effect of habitat loss on networks, are also 
important to conservation practice. Additionally, it has been 
suggested that the extinction of phylogenetically related species 
can lead to cascading coextinction events (Rezende et al, 
2007a). Understanding this potential effect with the assistance 
of network theory could help conservationists make decisions 
on species priority. 

Unfortunately, despite the theoretical advances in networks 
and their potential use in conservation, their actual 
implementation in conservation management is rare. There 
appears to be a communication failure among scientists, 
practitioners, and government officials that requires the 
assistance of all parties to resolve (Heleno et al, 2014). 
MUTUALISTIC NETWORK FROM AN 
EVOLUTIONARY PERSPECTIVE 


STUDIES 


As discussed above, studies on network structural properties 
and their implication for ecology have been extensively 
investigated in recent years. More attention has been spent on 
investigating mutualistic networks from the ecological 
perspective than the evolutionary one. However, species are 
not independent entities but rather related to each other through 
common evolutionary histories. Phylogenetic constraints may 
influence mutualistic interactions, imprinting a phylogenetic 
signal (i.e., the tendency of phylogenetically similar species to 
have similar phenotypic attributes; Bascompte & Jordano, 2007) 
on network structure (Ives & Godfray, 2006; Rezende et al, 
2007a, b). To understand plant-animal mutualistic networks, we 
must examine the phylogenetic signal of species’ positions (i.e., 
species centrality and their placement relative to modules; 
Olesen et al, 2007) in the network and observe which species 
form the network. 

In addition, studies that have attempted to link network 
literature to explore the evolution of networks (e.g., Guimaraes 
et al, 2011) remain scarce. The synthesis of network studies 
and phylogeny may change our understanding of the 
coevolutionary process. For example, recent work that 
combined phylogeny and pollination data found coevolution to 
be an important driver of species diversification (Van der Niet & 
Johnson, 2012). Despite some limitations of using phylogeny to 
explain the processes that influence speciation, work of this 
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kind may nevertheless suggest novel directions towards 
understanding speciation, diversification, and biodiversity. 
Below, we discuss three topics in the evolution of mutualistic 
networks in which there has been recent research activity. 


Phylogenetic signal in mutualistic networks 

A revolutionary article on network evolution used phylogenetic 
methods for the first time to study how the interaction pattern 
was associated with phylogeny in mutualistic networks 
(Rezende et al, 2007a). By incorporating a large dataset of 36 
plant-pollinator and 23 plant-frugivore mutualistic networks, they 
found that the phylogenetic signal in species degree (number of 
other species with which a species interacts) could be detected 
in more than one third of the networks. Meanwhile, the actual 
identity of interaction partners had a phylogenetic component in 
about half the interaction networks. Simulated extinction events 
triggered cascading coextinction, wherein phylogenetically 
related species went extinct together. 

Having realized the importance of evolutionary history on 
network structure, more researchers now include phylogenetic 
studies in their work. For example, a recent study by 
Schleuning et al (2014) has further improved our understanding 
on network modularity. These researchers associated both 
weighted and binary network data with species traits and 
phylogenetic information to study how modularity is related with 
these ecological and evolutionary factors. The study followed 
the methodology of Olesen et al (2007), which identifies species 
connectedness within (z-score) and between (c-score) modules. 
The results showed that for both weighted and binary networks, 
no phylogenetic signal was detected in within-module degree 
(z), but significant phylogenetic signal was found in the c value: 
the tendency of a species to interact with species in other 
modules. Results also showed significant phylogenetic signal in 
species degree, concordant with the results of Rezende et al 
(2007a). 


How do mutualistic networks affect coevolution of their 
species? 

Guimaraes et al (2011) combined a model for trait evolution 
with data from 20 plant-animal mutualistic networks to explore 
coevolution. They found that both evolution and coevolution 
contributed to the increase in convergence (e.g., trait similarities 
emerge in response to similar selective pressures) and 
complementarity (e.g., degree of trait matching between 
interactive partners) among traits within the network, and that 
coevolution significantly sped up the rate of trait evolution within 
networks (Guimaraes et al, 2011). They also found that super- 
generalists facilitated trait evolution in mutualistic networks by 
greatly increasing complementarity (match of traits between 
partners) and convergence. Because super-generalists serve 
as connections between different modules, evolutionary and 
coevolutionary forces that affect them multiply throughout the 
whole network (Guimaraes et al, 2011). 

Thompson (2009) proposed a new hypothesis from an 
evolutionary perspective, in which plant-animal interaction 
networks act as vortexes, absorbing increasing numbers of 
species into the interaction network. As the web grows bigger, it 


has more capacity to hold diverse species. In this way, plant- 
animal interaction networks may promote the evolution of 
biodiversity. However, more evidence is needed to support or 
reject this hypothesis. 


How could mutualistic networks affect the rate of 
evolution of participating species? 

The above studies either combined data on mutualistic 
networks with the phylogeny of the interacting species to look at 
how evolutionary history influenced network properties or 
simulated how coevolution shaped the patterns of 
complementarity and convergence between species in 
networks. However, the question of the rate at which species 
with different specialization levels evolve and diverge (i.e., 
speciation) in mutualistic networks remains largely unexplored. 
As has been well studied, mutualistic networks have a wide 
range of species that differ in their connectedness to other 
species, with many species interacting with only a few species 
(“specialists”), and other species interacting with many other 
species (“generalists”). In the field of classical ecological 
evolution, many studies have compared the speciation rate 
between specialists and generalists (e.g., Colles et al, 2009; 
Fernandez & Vrba, 2005); however, do different kinds of 
species in networks differ in their speed of evolutionary 
divergence? 

We recently conducted a study along the lines of this 
question, asking if specialist and generalist frugivorous birds 
have different speeds of evolutionary divergence (Gu et al, 
2015). A recent time-calibrated phylogeny of birds that included 
all extant species (Jetz et al, 2012) allowed such a study, as 
every species had an estimated time of divergence from its 
sister species. Using 16 seed dispersal networks in the 
published literature, we found that specialists (defined as per 
Olesen et al, 2007) had significantly shorter divergence times 
than did generalists. This result is somewhat surprising as most 
birds are thought to be rather unspecialized and not dependent 
on specific fruiting trees, in contrast to the very strong and 
sometimes obligate relationships found between pollinators and 
plant species (Blüthgen et al, 2007). We therefore treat the 
result quite tentatively, noting some confounding variables, such 
as specialists being generally rare, which could be driving the 
result. Nevertheless, we think the question is important, and 
believe future research should investigate this in other taxa 
such as insects and plants. 


FUTURE DIRECTIONS AND CHALLENGES 


A general challenge for studies of plant-animal mutualistic 
networks is the paucity of data sampling. Field observations are 
time-consuming, and certain habitats (e.g., canopies) are 
difficult to access. Fundamentally, there is a lack of properly 
trained personnel to make the observations. General problems 
for flora and insect taxonomy lie in the shortage of funding, the 
low Impact Factor index of taxonomic journals, low scientific 
regard of taxonomic knowledge, and a lack of investment and 
training programs for a new generation of taxonomists (Dangles 
et al, 2009; Ma, 2014; Rafael et al, 2009). These problems 


greatly hamper the development of plant-animal interaction 
network studies, leading to misidentification of species. 

There are potential alternative methods of network 
construction, including DNA barcoding (Carcia-Robledo et al, 
2013; Heber et al, 2004). For example, researchers recently 
found that DNA barcoding gave comparable results to direct 
observations in a plant-herbivore system (Carcia-Robledo et al, 
2013). However, the authors only studied a small group of well- 
known taxa (with plants only from order Zingiberales and rolled- 
leaf beetles from only two genera, Cephaloleia and Chelobasis). 
They emphasized that DNA information is unavailable in public 
databases for many organisms that participate in plant-animal 
interaction networks. This is especially true for mutualistic 
networks, in which the organisms come from a wide range of 
taxa. 

Additionally, the size of the network from different studies 
varies greatly, with total species ranging from a few dozen to 
several hundreds. This range in network size could reflect 
biological reality or could be due to differences in sampling 
effort. If the differences are due to sampling effort, these can be 
accounted for during analysis, as recently shown by Schleuning 
et al (2014). To answer large ecological or evolutionary 
questions, data from many different networks needs to be 
gathered together and a general protocol needs to be 
established (Heleno et al, 2014). 

Another challenge of network studies is to find an effective 
way to quantify the network, which should reveal the real 
ecological processes behind the interactions. For example, for 
most existing pollination networks and seed dispersal networks, 
flower visitation and frugivory are recorded as proxies without 
evaluating the effectiveness of the ecological services, that is, 
how these services affect plant reproduction (Heleno et al, 
2014). By reconsidering how networks are quantified, we can 
better understand the real ecological and evolutionary 
mechanisms. In addition, as weighted data are more 
informative than binary data (Barrat et al, 2004; Schleuning et al, 
2014), we should encourage future studies to use weighted 
network data. 

Further development of plant-animal mutualistic network 
theory requires more complete and informative datasets, which 
allow the evaluation of multiple mechanisms simultaneously 
(Vazquez et al, 2009). According to the geographic mosaic 
theory (Thompson, 2009), coevolution varies in time and space. 
Thus, network data with explicit spatio-temporal information are 
needed to provide more reliable and explanatory interpretations 
on both the ecological and evolutionary study of mutualistic 
networks. To answer questions of how networks affect the 
evolution of participating species, both comprehensive 
ecological information and phylogenetic information of 
participating species are required. We believe that field studies, 
particularly those that use the same systematic technique to 
look at networks in multiple study sites, and hence better 
understand consistency and variation in the qualities of the 
networks, are just as vital to the development of the field as 
theoretical studies. 

New insights will undoubtedly arise as our knowledge about 
phylogeny increases. In pollination networks, for example, there 
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are extremely large numbers of diverse pollinators from Diptera, 
Hymenoptera, Lepidoptera and Coleoptera. Yet only a small 
proportion of these species have been well studied 
taxonomically and phylogenetically. Time calibrated phylogenies 
for species in these groups are extremely limited (to some 
species from genus Bombus, family Apidae), and few of these 
species have DNA sequences accessible in GenBank. New 
developments in obtaining further phylogenetic information and 
better techniques will surely open new doors to analyses, just 
as Jetz et al (2012) did for our study (Gu et al, 2015). 

Despite all the challenges in mutualistic network studies, 
there are many opportunities in this field. It is difficult to predict 
the directions of future research to come, after all who could 
have foreseen the network literature blossoming as it has over 
the last 30 years, but certainly phylogenetic and bar-coding 
techniques may ignite new possibilities. From our evolutionary 
perspective, it is clear that there are many questions yet 
unanswered about how the network phenomenon has 
influenced the evolution of participating species. We also 
believe that conservation is a priority. However, despite the 
rapidly accumulating research in this field, there is little 
consensus on how environmental managers should incorporate 
networks into their planning (Heleno et al, 2014). For example, 
should networks be used to evaluate habitat quality rather than 
individual species counts (Valiente-Banuet et al, 2014)? Are 
networks where invasive species are  super-generalists 
irreparable or do remedial actions exist (Aizen et al, 2008; 
Valdovinos et al, 2009)? Further progress on these issues will 
demonstrate that our growing knowledge about ecological 
networks can be applied to solve environmental problems. 
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